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Abstract 

This paper is concerned with the stability problem for the planar linear switched system x(t) = 
u(i)Aix(i) + (l—u(t))A2x(t), where the real matrices Ai,A% € R 2x2 are Hurwitz and u(-) : [0, oo[— > {0, 1} is 
a measurable function. We give coordinate-invariant necessary and sufficient conditions on A\ and A2 under 
which the system is asymptotically stable for arbitrary switching functions u(-). The new conditions unify 
those given in previous papers and are simpler to be verified since we are reduced to study 4 cases instead 
of 20. Most of the cases are analyzed in terms of the function T(A\, A2) = ^(ti(Ai)tr(A2) — ti(AiA2)). 

Keywords: planar switched systems, asymptotic stability, quadratic Lyapunov functions 

1 Introduction 

Let A\ and A2 be two 2x2 real Hurwitz matrices. In this paper we are concerned with the problem of finding 
necessary and sufficient conditions on A\ and A2 under which the switched system 

x(t)=u(t)A x x(t) + (l-u(t))A 2 x(t), x = {xi,x 2 ) G K 2 , (1) 

is globally asymptotically stable, uniformly with respect to measurable switching functions u(-) : [0, oo[— > {0, 1} 
(GUAS for short, see Definition [T] below) . 

This problem has been studied in [3] in the case in which both A\ and A 2 are diagonalizable in C (diagonaliz- 
able case in the following) and in [T] in the case in which at least one among A\ and A 2 is not (nondiagonalizable 
case in the following). (See also [7] as well as the related work [5J.) 

In both cases the stability conditions are given in terms of coordinate-invariant parameters. Unfortunately 
the parameters used in the diagonalizable case become singular in the nondiagonalizable one and therefore the 
two cases were studied separately. 

The purpose of this note is to unify and simplify these conditions, reformulating them in terms of new 
invariants that permit to treat all cases at the same time. 

We have reduced the cases to be studied from 20 (14 in the diagonalizable cas^l and 6 in the nondiagonal- 
izable one) to the following 4 cases (see Theorem [T]). 



1 The first two authors were supported by a FABER grant of Universite de Bourgogne 

2 The stability conditions given in [3] were not correct in the case called RC.2.2.B. See [7] for the correction 
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SI: the first one corresponds to the case in which there exists a common quadratic Lyapunov function. The 
condition of SI is indeed equivalent to the condition given in |10j but is simpler to check. Recall however 
that the existence of a common quadratic Lyapunov function is only a sufficient condition for GUAS 
(i.e. there exist GUAS systems not admitting a quadratic Lyapunov function). See [11 [7] for details. 

S2: the second one corresponds to the situation in which there exists v e (0, 1) such that vA\ + (1 — v)A 2 has 
a positive real eigenvalue. In this case the system is unbounded since it is possible to build a trajectory 
going to infinity approximating the (non admissible) trajectory corresponding to u(t) = v and having the 
direction of the unstable eigenvector of vAi + (1 — v)A2. 

S3: in the third case there exists a nonstrict common quadratic Lyapunov function. The system is only 
uniformly stable, but not GUAS, since there exists a trajectory not tending to the origin when t goes to 
infinity. 

S4: in the fourth case the stability analysis of the system reduces to the study of a single trajectory called 
worst trajectory. If this trajectory tends to the origin then the system is GUAS (in this case there exists 
a polynomial Lyapunov function, but not a quadratic one). If it is periodic then the system is uniformly 
stable but not GUAS. If it is unbounded then the system is unbounded. 

For a discussion of various issues related to stability of switched systems, we refer the reader to [4j[5]. 

The paper is organized as follows. In Section 11.11 we recall the fundamental notions of stability and the 
different types of Lyapunov functions used in the paper. Section [2] contains our main result. In Section [3. II we 
define the normal forms that are needed in the proof. In Section 13.21 we give the details of the proof. 

1.1 Notions of stability 

Let us recall some classical notions of stability which will be used in the following. 

Definition 1 For 6 > let B$ be the unit ball of radius d, centered in the origin. Denote by hi the set of 
measurable functions defined on [0, oo[ and taking values on {0, 1}. Given xq € R 2 , we denote by 7 Xo ,u(-)(") the 
trajectory of ([I]) based in xq and corresponding to the control u(-). We say that the system ([I]) is 

• unbounded at the origin if there exist Xq £ M. 2 and u(-) £ U such that 7 KQjtl (.)(t) goes to infinity as t 
goes to infinity; 

• uniformly stable at the origin if for every e > there exists S > such that 7 X0 . «(.)(£) £ B s for every 
t > 0, for every it(-) £ U and every xq £ Bs; 

• globally uniformly asymptotically stable at the origin (GUAS, for short) if it is uniformly stable 
at the origin and globally uniformly attractive, i.e., for every 61,62 > 0, there exists T > such that 
7x ,«(-)(*) e for every t > T , for every u(-) £ U and every xq £ B$ 2 . 

Remark 1 The stability properties of the system {1} do not change if we allow measurable switching functions 
taking values in [0, 1] instead of {0, 1} (see for instance [7]). More precisely the system ([I]) with u(-) : [0,oo[— > 
{0, 1} is GUAS (rcsp. uniformly stable, resp. unbounded) if and only the system (Q]) with u(-) : [0, oo[— > [0, 1] 
is. In the following we name convexified system the switched system with u(-) taking values in [0, 1]. 

Since the stability properties of the system |T]) do not depend on the parametrization of the integral curves 
of Aix and A^x, we have the following. 

Lemma 1 If the switched system x = u(t)A\x + (1 — u{tf)Aix, it(-) : [0, oo[— > {0, 1}, has one of the stability 
properties given in Definition^ then the same stability property holds for the system x = u{t)(A\ / a\)x + (1 — 
u(t))(A 2 /a2)x, for every «i,a 2 > 0. 

Definition 2 A common Lyapunov function ( LF for short) V : IR 2 — > M + for a switched system of the form 
^ is a continuous function such that V(-) is positive definite (i.e. V(x) > 0, Vx ^ 0, V(0) = 0) and V(-) is 
strictly decreasing along nonconstant trajectories. 
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A positive definite continuous function V : M 2 — > K + is said to be a nonstrict common Lyapunov function 
if V(-) is nonincreasing along nonconstant trajectories. 

A common quadratic Lyapunov function (quadratic LF for short) is a function of the form V(x) = x T Px 
where P is a positive definite symmetric matrix and the matrices A[P + PA\ and A\P + PA2 are negative 
definite. 

We recall that, for systems of type |T]), the existence of a LF is equivalent to GUASjfl (see for instance [J]). 
Moreover the existence of a nonstrict LF guarantees the uniform stability of (Q} . 

2 Stability conditions for two-dimensional bilinear switched sys- 
tems 

We start this section by defining the notations and the objects that will be used to state our stability result. 
In the following the word invariant will indicate any object which is invariant with respect to coordinate 
transformations. As usual, we denote by det(X) and tr(X) the determinant and the trace of a matrix X. If 
X E R 2x2 the discriminant is defined as 

S x = tr(X) 2 - 4 det(X) . 
Given a pair of matrices X, Y we define the following object: 

T(X,Y) :=i(tr(X)tr(r)-tr(XF)). 

By means of these invariants we can define the following invariants associated to ([T]): 



Remark 2 Let us define 



tr(^i) 



V\Sa~\ 
2 

2nr 2 



if 5 Al ± , 5 A2 ^ , 
if 6 Al 5 A2 = but 5 Aj ^ , 
if 5a 1 = Sa 2 = , 
(tr(^M 2 ) - -tr(Ai)tr(A a )) , 



K := 



tr(Ai)tr(A 2 ) 
A := 4(T(A 1 ,A 2 ) 2 - T{A ll A 1 )T{A 2l A 2 )) . 

2r(Al ' A2) + ^ :e ^^, where, for i = l,2, 



2 v /det(Ai)det(A 2 



w — arctan 
arctanh 



tr(A 1 )tr(A 2 )(fcr,+r 3 _ I ) 
2riT 2 \/A 

2t!T 2 \/A 

tr( J 4 1 )tr(A 2 )(fer i -r 3 _ ! ) 
2-/A 



(tr(AiA 2 )-tr(A 1 )tr(A 2 )/2)r I 



if S Al < 
if S Ai > 
if 5 Ai = 



sign(x) := 



-1 if x > 
if x = 
-1 if x < 



Notice that, for every matrix X e R 2x2 , one has T(X,X) = det(X). Also, since the Killing form of R 2x2 
is defined as B{X,Y) = 4tr(XF) - 2tr(X)tr(F) one has sign(k) = sign{B{A 1 ,A 2 )). Finally, notice that if 
A\,A 2 are Hurwitz then r, < for i = 1,2 and signY{A\, A 2 ) — sign(r\T 2 — k). 



3 In [2 [3 \E\ Hj] it is actually shown that the GUAS property is equivalent to the existence of a polynomial LF 
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2.1 Statement of the results 

In this section we state our main result which characterizes conpletely the stability properties of two-dimensional 
bilinear switched systems. Our necessary and sufficient conditions apply both to the non-degenarate cases 
studied in [3] and to the degenarate ones studied in pQ. 

Theorem 1 We have the following stability conditions for the system ([1]) 

SI If T(Ai,A 2 ) > - v /dct(Ai) det(A 2 ), and tv(AiA 2 ) > -2 v /det(Ai) det(A 2 ) then the system admits a 



If - v /det(Ai)det(A 2 ) < T(A 1 ,A 2 ) < y/det(Ai) det{A 2 ) then the condition tr( A x A 2 ) > -2 A /det(Ai) det(A 2 ) 
is automatically satisfied. As a consequence the system admits a quadratic LF. 



S4 If T(A 1 ,A 2 ) > y /dct(A 1 )dct(A 2 ), and tr(AiA 2 ) < -2 v /dct(A i ) dct(A 2 ) then the system is GUAS, 
uniformly stable (but not GUAS ) or unbounded respectively if 

n<i,n = i,n>i. 

The following corollary will be derived from item SI of the previous theorem. 
Corollary 1 If det([Ai, A 2 ]) > then the system admits a quadratic LF. 

Remark 3 In the diagonalizable case Sa^a^ 7^ the parameters n, t 2 , and k are invariant under the transfor- 
mation (Ai,A 2 ) — > (Ai/ai, A 2 /a 2 ), for every ai,a2 > 0. This is no more true in the nondiagonalizablc case. 
Notice however that in any case the stability conditions of Theorem [T] do not depend on coordinate transfor- 
mations or on rcscalings of the type (A\,A 2 ) — > {A\/a\, A 2 /a 2 ). This is true in particular for the function 1Z. 



3 Proof of the main results 
3.1 Normal forms 

The aim of this section is to reduce all the possible choices of the matrices Ay,A 2 to suitable normal forms, 
obtained up to coordinates transformations and rescaling of the matrices (see Lemma [T] and Remark [3] above), 
and depending directly on the coordinate invariant parameters introduced above. The normal forms used here 
describe all the possible situations for two-dimensional bilinear switched systems, covering at the same time 
the diagonalizable case studied in [3] and the nondiagonalizable one studied in p] . They will play a key role in 
the proof of our results. 

Lemma 2 We have the following cases depending on the rank of [A\,A 2 ]: 

1. If det([.Ax, A 2 \) 7^ 0, up to a linear change of coordinates and a renormalization according to Lemma\I\ 
we can assume the following. 



quadratic LF. 



S2 If T(Ai,A 2 ) < —y/dct(Ai) dct(A 2 ), then the system is unbounded. 



S3 If T(A\,A 2 ) = — -y/det(Ai) dct(A 2 ), then the system is uniformly stable but not GUAS, 




(2) 




= 2k 




(3) 
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(b) 7/dct([Ai, A 2 ]) > then 5a, > for i = 1, 2, k e (-1, 1) and A 2 has the form 



. ( t 2 + ^1~ k 2 k A ... 

2. Under the hypothesis that rank{[A\, A 2 \) = 1, it is always possible, up to exchanging A\ and A 2 , to find 
a linear change of coordinates which diagonalizes A\ and renders A 2 upper triangular. 

3. If [Ai,A 2 ] = then it must be signlJiAf) = sign(SA 2 )- If °~Ai < for i = 1,2 or SA t > for £ = 1,2 
then, up to a linear change of coordinates and a renormalization according to Lemma Ql A\, A 2 assume 
the form @ and ([5]). respectively, with F = k = ±1. If Sax = o~A 2 = then A\,A 2 can be put in upper 
triangular form with the elements of Ai equal to Ti on the diagonal, for i = 1,2. 

Proof of Lemma [2j For simplicity we will prove the lemma just in the case det(L4i, A 2 ]) ^ 0, the other 
case being analogous. Note that Lemma [5] was proven in [3] when 5a x < 0, Sa 2 < an d in [T] in the case 
Sa 1 Sa 2 = 0. Therefore, we can assume either S Al > or 5a 2 > 0. First consider the case 5 Al > 0. In this case 
we can find a system of coordinates such that 

M = ( Q 1 M , A 2 = ( a c a,b,c,deR. (5) 

Without any loss of generality we can assume that Ai < X 2 . The discriminant of A\ is (A2 — A1) 2 and the 
discriminant of A 2 is 5a 2 = (a — d) 2 + Abe, which can be positive or negative. We have 

1^1= U\> 4( V 2) )- w 



Therefore detL4i,A 2 ] = bc(Xi - A 2 ) 2 and S\ A . A „i = -46c(Ai - A 2 ) 2 and k = -4=. If detL4i,A 2 ] < 



-\ 2 j aim u[ Au A 2 ] = -^H^l ~ ™ = " d ~ 

consider the linear transformation 

T 



— — \ 

v i r V ]■ 

which diagonalizes [A\, A 2 \. Then a straightforward computation shows that 

- 7 ==T-M 1 T=f 7 M and -^T^T = ( T l si 9<^)/F 

where -F satisfies the equation i* 1 + sign{5 a 2 ) / F = 2fc, and moreover we can assume \F\ > 1 up to eventually 
exchange the reference coordinates. If 6 Al < then it must be 5a 2 > and, exchanging the roles of A\ and 
A 2 , we can repeat the previous procedure obtaining 

A (n sign(S Al )/F \ ^ _ f 1 



F n / ' \ sign(5 A2 ) t 2 

Then the required normal forms are obtained by exchanging the coordinates and by a dilation along one of the 
coordinate axis. 

Consider now the case det [Ai, A 2 ] = 6c(Ai — A 2 ) 2 > 0. We have 

bc> => S A2 = (a - d) 2 + Abo and Ifel = , * Q = ^ <] 

y/{a - d) 2 + Abe 

In this case L4i,^4 2 ] is no more diagonalizablc. Using the transformation 

l^Jbc -lVb~c 



U = 
we get 



1 1 



2 -.U- 1 A 1 U=( T f M and -JL=U^A,U ( T2 + ^ 



which concludes the proof of the lemma. 
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3.2 Proof of Theorem CD and Corollary [T] 

To prove our main result we will assume, from now on, that A\,A 2 are under the normal forms given by 
Lemma [2j The following lemma, which can be proved by direct computation, will be used to take advantage 
of the conditions of [10] which describe the systems admitting a quadratic LF. 

Lemma 3 For any a £ [0, 1], we define 

4>{a) := det (aA ± + (1 - a)A 2 ) , ip{a) := dct (aA 1 + (1 - er)^ 1 ) . 

We have 

<P{ct) = cr 2 det A x +2cr(l -a)T{A 1 ,A 2 ) + (1 -cr) 2 det A 2 (7) 

and 

tP(a) = -r^-r-{(J 2 det Ai det A 2 + a(l - o-)tx{A 1 A 2 ) + (1 - a) 2 ) . (8) 



Proof of SI. Recall that the main result in [10] claims that the system {1} admits a quadratic LF if 
and only if <fi(a) > and ip{cr) > for every a e [0,1]. Notice from Lemma |3] that <p(cr) > if and only 
if either T(A\ 1 A 2 ) > or the discriminant A of JJJ) is negative. Analogously ip(o~) > if and only if either 
tr(AiA 2 ) > or the discriminant tr(AiA 2 ) 2 — 4det(Ai) det(A 2 ) of (jHJ) is negative. It is therefore clear that 
the cases considered in SI are those satisfying the conditions of [10]. The last statement of SI comes from the 
following series of inequalities 

1 1 



Vdet(Ai)det(A 2 ) > -(tr(Ai)tr(A 2 ) - tr(A x A 2 )) > --tr(AiA 2 ) . 

This concludes the proof of SI. 

Proof of Corollary [TJ To prove Corollary [T] in the case det([A 1; A 2 ]) > we use the point l.(b) of 
Lemma [5] In particular we have that T(Ai,A 2 ) = t\t 2 — k > and tr(AiA 2 ) = 2{t\t 2 + k) > so that the 
conditions of SI are satisfied. In the case det([Ai, A 2 ]) = the result was already known (see for instance [5]), 
and it can be easily proved by using the normal forms defined in Lemma [2] 

In what follows we will always assume det([Ai, A 2 ]) < 0. 

Proof of S2 and S3. Assume that T(Ai, A 2 ) < —y/det(A±) det(A 2 ). Then a straightforward computation 
shows that the minimum of </>(cr) is given by 

det A 2 - T(A U A 2 ) -A 
aa = dctA 1+ dctA 2 -2T(A u A 2 ) G ( °' 1} md ^ = 4(det A, + det A 2 - T(A U A 2 )) ~ °" 

In particular in the case described by S2 we have (j>(o~o) < and therefore the matrix ctqAi + (1 — <to)A 2 has a 
positive real eigenvalue, so that the system is unbounded (see Remark [T]). 

Similarly when A = and T(Ai,A 2 ) < we have <j>(o~o) — so that the system cannot be GUAS. In this 
case to prove that the system is uniformly stable it is possible to show that the system admits the following 
non strict quadratic LF: 



(sign(S Al )sign(S A2 ) - F 2 ) 2 
AF 2 (T Al F-T A2 sign(S Al )) 2 ' U ' 2 



v(x) = v( Xl ,x 2 ) = x 2 + ^ 2 



Proof of S4. First observe that, under the conditions of S4, we have F < — 1 and k < since, when A\ 
and A 2 are in normal form, tr(AiA 2 ) = F + sl ^ s p 5 -^) + 2l - lT2 = 2(k + nr 2 ) < and \F\ > 1. 
To prove S4 we introduce the set of points where the vector fields A-^x and A 2 x are parallel: 

Z = {x G R 2 : Q(x) = 0} , 

where Q(x) := det{A\x 1 A 2 x) . The discriminant of the quadratic function Q(x) coincides with A. Since A > 
then Z consists on a pair of noncoinciding straight lines passing through the origin. Take a point ieZ \ {0}. 
We say that Z is direct (respectively, inverse) in x if Aix and A 2 x have the same (respectively, opposite) versus. 
We have the following lemma. 
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Lemma 4 If Z is direct (resp. inverse) in xq £ Z \ {0} then Z is direct (resp. inverse) in every point of 
Z \ {0}. Moreover in the case S4 we have that Z is always direct. 

Proof of Lemma [4l Let Z = D\ U D 2 where D\, D 2 are straight lines passing through the origin. Let us 

observe that, if x £ Di 

3a,eK s.t. A 2 A^ 1 Axx = A 2 x = aiAtx , 

i.e. at is an eigenvalue of A 2 A~{ X and A\X belongs to the eigenspace associated to it. So a\a 2 = det(A 2 A^ 1 ) = 
dct(Ai) > ^ which implies that sign(ai) = sign(a 2 ) i.e. Z is either direct in every point or inverse in every 
point. On the other hand it is easy to verify that A^ 1 = de ^ A (2rild — A\), where Id denotes the identity 
matrix, which, in the case S4, implies 

a 1+ a 2 = tr(A 2 A^) = tr(^(2 n A 2 - A 2 A l} ) = > 

so that Z is direct. □ 

Let rrii be the slope of Di, for i = 1,2. Then, if vi is a vector spanning Di, the orientation of the vector 
A\Vi with respect to the radial direction is determined by the quantity 

sign(det (AiVi,Vi)) = sign(m 2 — sign(8 Ai )\ 

Similarly, the orientation of the vector A 2 Vi with respect to the radial direction is given by 

sign(det (A 2 Vi, Vi)) = sign ( miS ^9 n (^M) — J = sign(F 2 — mfsign{o~A 2 ))- 



Lemma 5 If Z is direct, i.e. if Y(A\, A 2 ) > 0, it must be 

sign(mf — signlS^)) = sign(F 2 — m 2 sign(SA 2 )) = +1 



Proof of Lemma [5j Since T(A±, A 2 ) > by the previous equalities we get that e := sign(m 2 — sign((5 J 4 1 )) = 
sign(F 2 — m 2 sign(o~A 2 )). If £ = we are in the conditions of SI, since [A\, A 2 ] =0. If e = — 1 then it must be 
szgn((5 J 4 1 ) = sign(dA 2 ) = 1- In this case we have F 2 < mf < 1 which is impossible since \F\ > 1. □ 

As a consequence the vectors AiX point in the clockwise sense for every x £ Z. This property allows to 
define the main tool for checking the stability of |T]) under the conditions of S4. 

Definition 3 Assume that we are in the conditions of S4 and under the normal forms of Lemma [H Fix 
Xg £ M 2 \ {0}. The worst-trajectory j Xo is the trajectory of ([l}, based at Xq, and having the following property. 
At each time t, 7 Xo (t) forms the smallest angle in clockwise sense with the exiting radial direction. 

Figure [1] expresses graphically the meaning of the previous definition. 

It is clear that the worst trajectory always rotates clockwise around the origin when <5^ i < for some 
i £ {0, 1}. If &Ai > for i = 1, 2 then the eigenvectors of A\ are (1, 1) T and (1, — 1) T , while the eigenvectors 
of A 2 are (1, F) T and (1, —F) T . In this case it is easy to check that mim 2 < and therefore, from Lemma[Sl 
without loss of generality we can assume 

F < m 2 < — 1< 1 < mi < —F. 

As a consequence D\ and D 2 divide the space into four connected components, each one intersecting the 
eigenspace of exactly one among A\ and A 2 . This implies that also in this case the worst trajectory rotates 
clockwise around the origin. This trajectory is the concatenation of integral curves of A 2 x from points of D\ 
to points of D 2 and integral curves of A\x from points of D 2 to points of D\ (see Figure [2]). 
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"worst trajectory": forms the smallest angle 
with the exiting radial direction 




Figure 1: The worst trajectory 



x = A,x'' 


m = " F m = m, 




~7^^ / / m=l 




/ ,'/ \ x = A, x 




\ \/ \m = -l 


k = A, x^ 


m = F 2 



Figure 2: Construction of the worst trajectory in the case 5a ± > 0, Sa 2 > 
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GUAS 



uniform stability 



instability 




7?<1 



72 = 1 



e 



(a) 



(b) 



(c) 



Figure 3: The worst trajectory: meaning of the function 72 



As explained in the previous papers [TJ [3J the behaviour of the worst trajectory is sufficient to derive the 
stability properties of JT]). Let us analyse the worst trajectory 7 Ko (-) where x € A- Assume that T > is 
such that x\ = 7z (T) is the first intersection point between the worst trajectory and Di. The worst trajectory 
tends to the origin as time goes to infinity if and only if 72 := |xi|/|a;o| < 1, and in this case the system is 
GUAS (see Figure [3] (a)). It is periodic if and only if 72. = 1, and in this case the system is uniformly stable 
but not GUAS (see Figure |3J (b)). It blows up if and only if 72 > 1, and in this case the system is unbounded 
(see Figure |3J (c)). 

The computation of 72 was done in details in [TJ [3J . The formula which is given in Section [5] is a simpler 
reformulation of the ones previously obtained, in terms of our invariants. This concludes the proof of S4. 
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